Phase boundary and finite temperature crossovers of the quantum Ising model in two dimensions 
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We revisit the two-dimensional quantum Ising model by computing renormalization group flows close to 
its quantum critical point. The low but finite temperature regime in the vicinity of the quantum critical point 
is squashed between two distinct non-Gaussian fixed points: the classical fixed point dominated by thermal 
fluctuations and the quantum critical fixed point dominated by zero-point quantum fluctuations. Truncating 
an exact flow equation for the effective action we derive a set of renormalization group equations and analyze 
how the interplay of quantum and thermal fluctuations, both non-Gaussian in nature, influences the shape of 
the phase boundary and the region in the phase diagram where critical fluctuations occur. The solution of the 
flow equations makes this interplay transparent: we detect finite temperature crossovers by computing critical 
exponents and we confirm that the power law describing the finite temperature phase boundary as a function of 
control parameter is given by the correlation length exponent at zero temperature as predicted in an e-expansion 
with e = 1 by Sachdev, Phys. Rev. B 55, 142 (1997). 
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The quantum Ising model serves as a prime textbook ex- 
ample to illustrate fundamental aspects of quantum phase 
transitions ! 1 ' 2 ' 3 ' 4 ' 5 The quantum Ising Hamiltonian has the 
form, 



Hi 



Qi 



(i) 
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where J is a ferromagnetic exchange coupling, the sum (if) 
runs over pairs of nearest neighbor sites, and the quantum de- 
grees of freedom are represented by the operators erf 1 which 
reside on a site i of a hypercubic lattice in d dimensions 
and reduce to the Pauli matrices in the basis where cr is 
diagonal. 2 The parameter h is the external transverse magnetic 
field which induces quantum-mechanical tunneling events that 
flip the orientation of the Ising spins. The relevant parameter 
of Eq. (Q]l is the ratio 5 ~ J/h. For large 5 the ground state is 
ferromagnetically ordered and spontaneously breaks the dis- 
crete Z2 Ising symmetry while for smaller 6 the spins in the 
ground state remain disordered. The two phases are separated 
by a second order quantum phase transition at a critical <5 C . At 
finite temperature the formation of spin order is hindered and 
the 5 at which the order sets in is increased leading to a line 
of second order phase transitions T c (&j that terminates at the 

quantum critical point (QCP) T c (S C J = 0. Since the phase 
diagram of the quantum Ising model exhibits many generic 
features of physical systems in vicinity of their QCPs, it is 
important to understand it in detail. 

Various finite temperature properties of compounds mod- 
elled by the quantum Ising model were measured experimen- 
tally in three dimensions. 6 Theoretically, the corresponding 
phase diagrams were investigated by Sachdev within analyti- 
cal approaches. 27 These rely on the effective continuum field 
theory to which an expansion around the upper critical dimen- 
sion is applied. In two dimensions, the quantum Ising model 
describes a strongly coupled lattice system. Its ground state 
was recently analyzed numerically with new algorithms. 8,9 
The perturbative renormalization group (RG) approach by 



HertaiS and Millisii does not cover this case as the QCP is 
associated with a non-Gaussian fixed point therefore invali- 
dating the - in other cases successful - expansion around a 
Gaussian fixed points 

In this note, we extend our recent RG approach 12 to QCPs 
associated with non-Gaussian fixed points. We present an 
analysis addressing the quantum Ising model in two spatial 
dimensions near the QCP with flow equations derived within 
the functional renormalization group framework. 13 This set 
of coupled differential equations is valid at zero and finite 
temperature and is derived from a truncation of the exact 
functional flow equation for the scale-dependent effective ac- 
tion r[<f>], with (p a scalar-valued bosonic field obtained from 
coarse-graining Ising spins over a neighborhood of their lat- 
tice sites. The solution of this flow equation as a function of 
the continuous cutoff scale A yields the renormalized effective 
action from which all physical properties can be extracted. Al- 
ready in simple truncations, this framework yields the critical 
properties of 0(7V)-symmetric field theories below the upper 
critical dimension, 13 including the Ising case (9(1).— 

The scale-dependent action parametrizing the continuum 
field theory for the low energy physics of the quantum Ising 
model, 2 7 which we will apply in this note is given by 



P z(w 2 + p 2 )<f>- p + U[<f>] . (2) 



Here p = (o>„,p) and a>„ = 2nnT with n integer are bosonic 
Matsubara frequencies, Z is a A-dependent renormalization 
factor multiplying the momentum dependence of the propa- 
gator, and U[(p] is the effective potential specified below. For 
simplicity, the renormalization factors corresponding to the 
p 2 and col terms are taken to be equal here. As we checked 
by an explicit calculation, considering two different Z-factors 
for these two terms has no impact on the results of this pa- 
per. The action is regularized in the ultraviolet by restricting 
momenta to |p| < Auv- By virtue of the quadratic frequency 
dependence, the dynamical exponent z is equal to unity. The 
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effective dimensionality at zero temperature D = d + z = 3is 
below the upper critical dimension £) + = 4. We approach the 
phase boundary and the QCP from the symmetry-broken re- 
gion in the phase diagram and we therefore assume a potential 
U[<p] with a minimum at a non-zero order parameter </>q: 



r l/T r 

= J dr J d*jc 



4! 3! 2! 



(3) 



The 



where <p and (/>' are functions of x and r with (f> = <po 
parameter 5 - ^p- is related to the ratio of transverse field to 
exchange coupling 6 ~ J/h in the original Hamiltonian, Eq. 
([T), and controls the distance from criticality. The three-point 
vertex y3ud generates an anomalous dimension of the order 
parameter field already at one-loop level. 

The flow equations are obtained along the lines of Ref. 
The recipe is the following: after endowing the propagator 
with a suitably chosen cutoff function (also given in Ref. [l2l) 
that implements the A-dependence and regularizes the in- 
frared singularity of the massless propagator at criticality, one 
executes a cutoff-derivative on the analytic expressions cor- 
responding to all one-loop one-particle irreducible Feynman 
diagrams for the parameters u, <po. The flow of Z is obtained 
by taking second derivative of the equation describing the flow 
of the propagator with respect to momentum. After utilization 
of the variables 



b\Zd 



2KjTA d ~ 2 



u 2K d T 
dZ 2 A 4 - d 



(4) 



where is defined via j ^7 = Kd f d\k\ \k\ d as well as 



the anomalous dimension 77 and rescaled temperature T 



dlogZ 
d\ogA 



T = 



2nT 
~~A 



(5) 



we can write the flow equations as 



dli 



a! log A 



(d - 4 + 277) u + 



3ft 2 



dp 



a! log A 



+ -((«f) 2 + i + f) 

(2-d-T])p + 



1 y 1 



(6) 



The anomalous dimension is determined by 
1 2 



77 = 2ft 2 p 



+ 



2 I 



(l + ?f) (rf + 2)(l + 2 f) 

1 2 



"=1 {[nff + 1 + ?f ) 4 (d + 2) ((nT f + 1 + 



(7) 



The Matsubara sums can be performed analytically yielding 
hyperbolic functions but the expressions do not deliver any 
additional insights here. There are three distinct contributions 
in the flow equations Eq. (|6]): the Gaussian terms linear in p, 
ft, the classical terms corresponding to the zeroth Matsubara 
frequency, and the quantum terms summing over the non-zero 
Matsubara frequencies that become crucial at low tempera- 
tures. The zero temperature limit of Eq. (|6]l results in mod- 
ified variables p and ft, which then do not depend on T and 
are rescaled by different powers of A, but is fully accessible 
within our frameworkji 2 ^ A strength of the present approach 
lies in the ability to integrate out classical and quantum fluc- 
tuations - both non-Gaussian in nature, including the anoma- 
lous dimension rj of Eq. 

The numerical solution of the coupled flow equations (|6] 
13 at T = gives access to the critical behavior at the QCP 
while for T > the classical critical behavior along the phase 
boundary is obtained. As initial conditions for the upper cut- 
off, the interaction coupling, and the momentum renormaliza- 
tion factor we set Auv = 1, u = 1, and Z = 1. By choos- 
ing the initial condition of p or 5, we can tune to the critical 
state characterized by 6 vanishing at the end of the flow, that 
is, for A — > 0M- In Fig. Q] the zero temperature flow at the 
QCP is juxtaposed with finite temperature flows at the phase 
boundary. Both, the finite-T and the T = theory find a de- 
scription in terms of two distinct non-Gaussian fixed points 
indicated by two distinct scaling plateaus with finite ft and 77. 
The anomalous dimension directly at the QCP (black dots in 
Fig. [T](c)) comes out as 77QCP ~ 0.1 to be compared with 
the accurate value of the classical Ising universality class in 
3D:— 77 ss 0.04. The correlation length exponent at T — 
comes out as v = 0.6 in our calculation - coincidentally close 
to the currently accepted value v = 0.63. 18 We calculate v by 
first computing the susceptibility exponent y and then using 
the scaling relation y — v (2 - 77). 19 The susceptibility corre- 
sponds to the renormalized value of S^\ Q at the end of the 
flow. The exponent y is obtained as the slope of the graph of 
<5 A ^ versus the distance from the QCP (6 - 5 c )a=a uv m dou- 
ble logarithmic coordinates. 

At the finite temperature fixed point, 77 * 0.5, to be com- 
pared with the exact value from the Onsager solution of the 
classical Ising model in two dimension 77 = 1/4. 19 More elab- 
orate truncations lead to improved accuracy in the critical 
exponentS i 14 i 15 ' 16 

From the finite-temperature flows in Fig. [T](a), we can de- 
duce the Ginzburg-scale, where ft starts to become sizable, to 



vary with temperature as Ac °c T t 



l/(4-d) 



with d — 2, thereby 



fitting the formula valid for d + z > 4. Juxtaposing the 
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FIG. 1 : (Color online) Flows for the quantum Ising model in d = 2 as 
a function of logarithmic cutoff-scale s = - log [A/A uv ] for various 
temperatures and at T = 0. We set Auv = 1. The infrared (ultravi- 
olet) is to the right (left) of the graphs. The values of the classical 
fixed point, attained by all finite-r flows, and the quantum critical 
fixed point, attained by the zero-temperature flow, are marked on the 
vertical axis, (a): Flows of the quartic self-interaction u. (b): Corre- 
sponding flows of the rescaled minimum of the effective potential p. 
(c): Flows of the anomalous dimension rj. 



Ginzburg scale with the quantum-to-classical crossover scale 
which follows from the definition of f in Eq. (|5): A c / ~ T l ^ z 
with z — 1 , we obtain 



.1/2 



> A r 



(8) 



which indicates that non-Gaussian fluctuations become im- 
portant at energy scales above the quantum-to-classical 
crossover. For the cases where the QCP is described by a 
Gaussian fixed point, 12 this relation is inverted Ac < A c ;, 
while Ac ~ A c / for d = 2, z = 2. 
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FIG. 2: (Color online) Emergence of three different regimes in the 
phase diagram as illustrated with three different values of the order 
parameter exponent /? upon approaching the QCP at very low tem- 
peratures. Here, we set T = 1.7 x 10 5 . 



Another important difference between the QCP being Gaus- 
sian or non-Gaussian is that in the latter case Ac does not 
vanish as T — > 0. In the (T, S)-plane there is a jump from 
the 2D-Ising fixed point at finite temperature to the 3D-Ising 
fixed point at T — with a finite Ap =0 as shown by the zero- 
temperature flow in Fig. Q](a). 

The phase diagram in the symmetry-broken phase at low 
but finite temperatures can be portioned into three regimes, 
characterized by exponents belonging to different universality 
classes as shown for the order parameter exponent /3 in Fig. 
[2] and for the anomalous dimension 77 as a function of scale 
in Fig. [3] The regimes and their relative size in the phase 
diagram can be detected by computing /3 when approach- 
ing the phase boundary from the symmetry-broken phase: 
<po ~ yfp ~ (S — 5cf. The results in double logarithmic co- 
ordinates, so that the slope corresponds to the exponent ft, are 
exhibited in Fig. [2] In the immediate vicinity of the phase 
boundary, /3 ~ 0.13, which comes out close to the exact value 
/3 = 1/8 = 0.125 from the Onsager solution of the classical 
2D-Ising model. 19 This reflects the fact that at any finite, even 
if small, temperature the asymptotic properties of the system 
are determined by strong classical, non-Gaussian fluctuations. 




FIG. 3: (Color online) Crossover behavior of the anomalous dimen- 
sion as a function of scale s = - log [A/A uv ] for T = 1 x 10~ 6 . Upon 
increasing T, this curve will continuously deform toward the shape 
of the blue dashed-dotted line in Fig.[TJ(c). 
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FIG. 4: (Color online) Schematic plot of the Ginzburg-line (dashed, 
blue), the crossover line (dotted, black) from 3£>-Ising to 2£>-Ising 
behavior, and the true r r -line (straight, red) for the quantum Ising 
model in d = 2. Units are arbitrary and will depend on the mi- 
croscopic details of the specific system under investigation. The 
symmetry-broken phase is the area to the right of the r c -line. 



This asymptotic behavior does not survive the limit T — > 0, as 
the classical fluctuations are replaced by zero-point quantum 
fluctuations. Further away from the phase boundary, in the 
center of Fig. [2j /3 » 0.29, which comes out close to the value 
of the 3Z)-Ising universality class f3 = Q.33M This regime per- 
sists at zero temperature and reflects the non-Gaussian char- 
acter of the QCP. Further away from the phase boundary, to 
the right of Fig. |2l mean-field behavior sets in with /3 ~ 0.5. 

The same crossover between the three regimes manifests it- 
self also in the scaling behavior of the propagator represented 
by the anomalous dimension in Fig. [3] In the high energy 
regime of the flow (0 < s < 2), 77 is close to zero reflecting 
mean-field behavior. At lower energies, between s ~ 5 and 
s ~ 12, the flow is governed by the quantum critical fixed 
point which belongs to the 3d-Ising universality class with a 
scaling plateau at 77 ~ 0.1. Asymptotically in the infrared 
(15 < s < 18), classical scaling of the 2d-Ising universality 
class sets in with 77 being attracted toward 77 ~ 0.5. Upon in- 
creasing T, the size of the classical 2d-Ising plateau increases 



and the size of the quantum 3d-Ising plateau decreases. 

By identifying the finite temperature phase boundary with 
the critical coordinates in the (S, T)-plane where the non- 
Gaussian scaling plateaus in Fig. [T]occur, 12 we have computed 
the critical temperature as a function of the control parameter 
and obtained the power law 

T c ~(S- S c f (9) 

with v = 0.6 the correlation length exponent at zero temper- 
ature as determined above. Sachdev obtained the same re- 
sult for T c expanding around the upper critical dimension in 
e - 3-d and then extrapolating the result to d — 2. 7 Although 
yielding the correct result, at least the formal justification of 
setting e = 1 without capturing the anomalous dimension and 
the associated non-Gaussian behavior in the strong-coupling 
region of the phase diagram is not immediately clear to us. 

In Fig. |U we schematically plot the r c -line, the crossover 
line from 3D-Ising to 2D-Ising behavior, and the Ginzburg 
temperature (at which the interaction u becomes sizable). Re- 
lying on the Ginzburg-line as a proxy for the phase boundary 
gives a different location of the QCP. 

In conclusion, we have extended our recent RG frame- 
work for quantum-critical systems with discrete symmetry- 
breaking 12 to systems where the QCP is associated with a 
non-Gaussian fixed point. As a stresstest, we performed an 
RG analysis of the low but finite temperature regime of the 
two-dimensional quantum Ising model where the system is 
squashed between the strong-coupling fixed points of (i) the 
QCP and (ii) the finite temperature phase boundary. 

Interesting avenues for future investigation include the as 
yet unresolved interplay of quantum and thermal fluctuations 
in the vicinity of a QCP where the finite temperature phase 
transitions are of the Kosterlitz-Thouless type. 
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